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Abstract 



The azimuthal asymmetry of J/ip suppression in non-central heavy- ion 
collisions is studied within a dynamic model of J/tp suppression in a de- 
confined partonic medium. Within this model, J ftp suppression in heavy-ion 
collisions is caused mainly by the initial state nuclear absorption and dissocia- 
tion via gluon-J/?/> scattering in deconfined partonic medium. Only the second 
mechanism gives arise to azimuthal asymmetry of the final J/ip production. 
We demonstrate that if there is an onset of suppression by quark-gluon plasma 
(QGP) in the NA50 data, it must be accompanied by the non-vanishing az- 
imuthal asymmetry. Using the same critical density above which the QGP 
effect enters, we predict the azimuthal asymmetric coefficient V2 as well as the 
survival probability for J/ijj at the RHIC energy. 

PACS number(s): 12.38.Mh; 24.85. +p; 25.75.-q 
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In the search for quark-gluon plasma (QGP), J /if) suppression has been proposed as 
one of the promising signals [Q of the deconfinement in high-energy heavy-ion collisions. 
Because of the color screening effect in a quark-gluon plasma, the linear confining potential 
in vacuum that binds two heavy quarks to form a quarkonium disappears so that it can 
be easily broken up causing suppression of the J /if} production. The problem in heavy- ion 
collisions is however complicated by other competing mechanisms such as initial nuclear 
absorption |2| and hadronic dissociation |J . While recent precision data from the NA50 || 
experiment at the CERN SPS energies clearly show anomalous suppression unexplained by 
the normal initial nuclear absorption, there are still much debates about the exact nature of 
the anomalous suppression PU, whether it is caused by the formation of QGP or dissociation 
by ordinary hadronic matter. 

We propose in this letter the study of azimuthal asymmetry of J /if) production [|7| as 
additional measurements to distinguish different competing mechanism of J / if> suppression. 
Since the initial state interactions such as nuclear absorption or nuclear shadowing of gluon 
distribution has no preference over the azimuthal direction they will not have any contribu- 
tion to the azimuthal anisotropy of the J / if) production. Only suppression by the final state 
interaction with the produced medium will cause significant azimuthal anisotropy in the 
final J / if) distribution in the transverse direction. If the centrality dependence of the J/ if) 
suppression additional to the initial nuclear absorption is caused by formation of QGP, it 
must be accompanied by a sudden onset of the azimuthal anisotropy. On the other hand, a 
hadronic absorption scenario would give a continuous centrality dependence of the azimuthal 
anisotropy. In this letter, we will study the centrality dependence of both the averaged J/ if) 
suppression factor and the azimuthal anisotropy with a model in which J /if) suppression is 
caused by initial nuclear absorption and final state dissociation by QGP above a critical 
density. Using parameters from fitting the NA50 data, we will also give predictions for J/ if) 
suppression and its azimuthal anisotropy at the RHIC energies. 

We follow the microdynamic approach of J/ if) suppression M , in which the J/ if) suppres- 
sion is caused by gluonic dissociation. The gluonic dissociation cross section in this model 
depends strongly on the gluon density and the averaged energy of the gluons in the medium. 
In a hadron gas, one finds the effective hadron- J/-?/> cross section to be extremely small 
because of the small gluon distribution inside a thermal hadron. Furthermore, average en- 
ergy carried by the gluon inside a thermal hadron is much smaller than the threshold energy 
of J /if) dissociation. This further suppresses hadron- J /if) dissociation cross section. We 
should emphasize that such an assumption is model dependent. Results of other theoretical 
calculations |9|-[T2[| of hadron- J /if) dissociation also vary widely, generally potential models 
giving large cross sections while effective interaction models giving much smaller values. For 
consistency we will assume a small hadron- J /if) cross section within the gluon dissociation 
model. Therefore, we can neglect the J /if) dissociation in a hadron gas which exist in the 
peripheral heavy-ion collisions and the late stage of central collisions. 

Different from a normal hadron gas, a deconfined partonic system contains large number 
of much harder gluons which can easily break up a J /if). This will be the only mechanism 
for final state suppression of J /if) in our paper. Such a gluon dissociation model can be 
considered as a micro-dynamic model of screening of the potential between heavy quarks 
in a quark-gluon plasma J|.[I^.[13|] . In the screened potential model, a J /if) is considered 
destroyed if it is produced in a plasma that is above the critical temperature. On the other 
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hand, the micro-dynamic model still gives a non-vanishing survival probability that depends 
on the temperature of propagation length. 

The perturbative calculations predict the gluon- J/ip dissociation cross section [ YA , FJ ] . 



where 

AT ^(M,^. 

3 3 3g 2 s m 2 Q 

Here # s is the strong coupling constant, tuq is charm quark mass, and q° is the gluon energy 
in the J/ip rest frame. To break up a J/ip, q° must be larger than the binding energy e . 
Using similar approach, we have calculated the dissociation cross section for P-wave states 
by gluons, 

where e x is the binding energy of the P-wave state. For Xc it is about 0.250GeV. Because 
of the small value of the binding energy, the validity of the perturbative calculations for the 
gluonic dissociation cross section of Xc might be questionable. So the above formula can only 
be considered more phenomenological. However, keeping this point in mind, we can see that 
the above expression still qualitatively reflects the fact that Xc states are easier to be broken 
up than J /if), because of the much lower energy threshold and the overall larger factor of 
4. Therefore, this perturbative calculation gives us a reasonable estimate and guides us to 
include its contribution for a more complete study of J / if) suppression in heavy ion collisions. 
In the following we will consider Xc contributing to about 40% of the initial J/ip production 
and use the above formula to estimate its suppression in a deconfined partonic system. 

In the rest frame of a deconfined parton gas, the momentum distribution of thermal 
gluons will depend on the effective temperature T with an approximate Bose-Einstein dis- 
tribution, f(k°;T) oc [exp(A;°/r) — l] -1 . The velocity averaged dissociation cross sections 
for charmonia is defined as, 

M(T,p) = Jd3kf{k0 . T) , (3) 

where v re i is the relative velocity between J/ip and a gluon. In this paper we are only 
interested in J/ip production in the midrapidity region, so these cross sections will depend 
on charmonia transverse momentum as well as the effective temperature T. 

With the velocity averaged dissociation cross sections, the survival probabilities of char- 
monia in the deconfined quark-gluon plasma will have the following form, 

gdeconf.^ ^ ^ = exp{ _ H ( Vrel(r)p ^ + tf T) r)0 ( p _ pc )dT}, (4) 

J TO 

where To is the formation time of the quark-gluon plasma, which will be set as r = lfm 
in the following calculations. The upper limit of the time integral Tf is determined by the 
function. Here we introduce the critical density p c above which the QGP dissociation 
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effects enters. We assume it is the same for both J/ip and Xc- Because of different binding 
energies, the effective cross sections will have different temperature dependences for J/ijj and 
Xc even if p c is reached, p is the local density depending on r, the initial production point 
r and the velocity v of the charmonium particles. The velocity v depends on the transverse 
momentum and the azimuthal angle <fi. The dissociation cross sections (v re ia) depend on 
the effective temperature, which will also depend on the local density p T (b,s). In the case 
of 1+1D Bjorken longitudinal expansion with the initial plasma density po = p( r = r o), 

p T (b,s) = p (b,s)(^r, (5) 

T 

where a = 1. Correspondingly, for the effective temperature, 

t t (M = 7UM(-) 1/3 - (6) 

T 

For simplification, we relate the local effective temperature with the local density in the 
following way, 

T T (b,s) = K(p T (b,s)) 1/ \ (7) 

where k = [tt 2 /16C(3)] 1/3 . 

The initial plasma density is related to the rapidity density of gluons, 

„ s 1 dN 9 . . 

where d 2 s is the transverse area of the overlapping region of two colliding nuclei. We will 
follow the two-component model |16j and include both the soft and hard contribution to the 
final hadron production. Assuming that the initial gluon density is proportional to the final 
hadron rapidity density, we have phenomenologically, 

dN 9 

(b, s) = c[f b n b (b, s) + f p n p (b, s)], (9) 



dyd 2 s 

where c is a constant and we will set c = 1 in the following calculations, and 

n p (b, s) = T A (s)[l - exp(-T B (b - s)a pp )] + T B (b -s)[l- exp(-T A (s)a pp )), (10) 
n b (b, s) = T A (s)T B (b - s)a pp . (11) 

Since mini-jet cross section at the SPS energy is very small, we will effectively only have the 
soft contribution which is proportional to the number of participant nucleons. At collider 
energies such as RHIC the contribution from the hard processes is more important. We will 
use fb = 0.34 and f p = 0.88 as determined by the PHENIX [0 experiment for Au + Au 
collisions at yfs = 130 GeV. We extrapolate these parameters to = 200GeV by just 
multiplying a factor of 1.14 found by PHOBOS pp . 

The final expression for the survival probability due to QGP suppression is, 

-* f r f dT 

S decon f(b,r,v)=exp{- / —(v rel a^ g )(f+vT,T)p Q (b,f+vr))Q(p - p c )}. (12) 

J to T 



4 



It is interesting to note that with a very large constant dissociation cross section the above 
formula will be equivalent to the model || by Blaizot et ai, where they assume that all 
of J/ip will be dissociated above some critical density. The detailed comparison of this 
approach and ours will be presented elsewhere. 

We can also include the transverse expansion effects on the local parton density following 



Ref. m, 

p(r, t) = — J ^Y^Poir- v t r), (13) 



where v t is the average velocity of the transverse expansion of the parton system. We will 
use v t = 0.4c for SPS and v t = 0.6c for RHIC in the following numerical calculations. 

Apart from the above discussed QGP suppression for charmonia states, there is also 
suppression associated with the initial state interaction, i.e., the nuclear absorption of so- 
called preresonance of cc pairs, 

S abs (b f) = ( X ~ exp(-(r afca T j4 (r)))(l - exp(-q afcs T B (|& - f|))) 

° 2 abs TA(r)T B (\b - r\) 

where a a bs is the absorption cross section of the preresonance with nucleons, for which we 
will set a^s = 5.8mb in this paper. 

By summing up these two contributions, we get the final J/ip survival probability as 

£ = J d*rT A (r)T B (\b - f\)S abs (b, f)S d ^jr, g) ^ 
f d*fT A (r)T B (\b - r\) 

From the above expression, we see that the nuclear absorption has no dependence on the 
azimuthal angle. This means that there is no contribution to azimuthal anisotropy from the 
initial nuclear absorption. On the other hand, the final state interaction or the dissociation 
by the deconfined parton gas indeed has azimuthal angular dependence because the parton 
density is azimuthally asymmetric. Therefore, any finite value of angular anisotropy for 
J/i/j production should come from the final state interaction. It will provide us important 
information about the early stage of the quark-gluon plasma. 

We will quantify the azimuthal anisotropy by the second Fourier coefficient v 2 of the 
azimuthal angle distribution of the final J / ip distribution, similarly to the proposed elliptic 
flow measurement pO| , pT| . It is defined as 



Jd4>S sur - (bjr cos 20 
W,Pt) = — — ; — t^— , 16 

where <fi is the azimuthal angle between J/ip transverse momentum and the impact 
parameter b. With this formula we can study both the px dependence and the centrality (b) 
dependence of v 2 . 

In the following we present the numerical results of the above approach. We first deter- 
mine the critical density p c by fitting the SPS data on J/ip suppression, and then predict 
the suppression and v 2 at RHIC. 

In Fig. 1, we show our results at the SPS energy. The upper plot is the survival probability 
as a function of transverse energy Et compared with the experimental data from NA50 [H . 
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From the fit we determined the critical density p c = 3.3/m 3 . The correlation between the 
impact parameter b and transverse energy Et p2 |, 



P(E T ,b) = =exp 



[Et - giVp 



27iq 2 aN p {b) I 2g 2 aiV P (&) J ' 

has been used in the calculations, where we set the parameters as q = 0.274GeV and 
a = 1.27 0. In the calculations, we also include the Et fluctuation effects as in [pU|p3"|l, 
which is important for the last few Et bins. From this figure, we can see that our approach 
can well describe the experimental data of NA50, and taking into account the transverse 
expansion improves the fit especially in the peripheral region. As expected, t>2 vanishes 
at very peripheral collisions, and becomes sizable when the anomalous suppression makes 
sense. For more central collisions, because of the symmetric geometry of the collisions, the 
azimuthal asymmetry v<i vanishes again. Since we assumed that only final state interactions 
produce finite value of t>2, its increase with Et is quite abrupt around the value when the 
critical density p c is reached. After taking into account the transverse expansion, t> 2 is a 
little higher than the case without transverse expansion. This is quite different from w 2 for 
jet quenching where transverse expansion is found to reduce v-i significantly. 



The results at RHIC are shown in Fig. 2. The survival probability and of J /if) are 
shown as functions of number of participants. Comparing with the results in Fig. 1, we find 
that the anomalous suppression enters already at peripheral collisions at RHIC, and the gap 
between the full suppression and the nuclear absorption alone is much larger than that at 
SPS. And i>2 is also much larger. 

It is interesting to note the different dependences of v-i on the the transverse expansion 
at the SPS and RHIC energies. At SPS, it enhances v-i a little bit while at RHIC it reduces. 
This is because at these two different energies the dominant suppression sources are different 
in the region where v-i is sizable. At SPS, a large part of J /if) suppression comes from Xc sup- 
pression, because the density is not so high and it is difficult to dissociate directly produced 
J/if). For this part, the dissociation cross section is very large, in the order of a few mb. So 
Xc would be dissociated almost totally when the local density is above the critical value p c , 
which is more similar to the Blaizot-like model. In this case, the transverse expansion will 
increase v-i a little bit because it increases the time duration when Xc is being dissociated. 
This increase of time duration even overcomes the decrease of geometrical asymmetry due to 
transverse expansion leading to an increased i^. The dominant suppression, however, comes 
from directly produced J /if) (about 60%) at the RHIC energies because of the much higher 
initial density. The transverse expansion accelerates the decrease of the initial density and 
reduces the initial geometric asymmetry. This leads to the decrease of t>2, very similar to 
the case of jet quenching in dense. Since the final total J /if) suppression is a mixture of 
direct suppression and suppression through Xc, this reduction is not as large as it is for jet 
quenching [ pjj] . 



In conclusion, we have studied the azimuthal asymmetry of J /if) suppression at both 
SPS and RHIC energies within a dynamic model of charmonia dissociation in a deconfined 
partonic system. We assume the hadronic absorption in our model is negligible so that the 
only final state J /if) absorption is due to gluon dissociation in deconfined medium above 
a critical parton density. With a critical density p c = 3.3/m -3 we can reproduce well 
the anomalous suppression found by the NA50 experiment at the SPS. We predicted the 
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azimuthal anisotropy t> 2 of the J/ip suppression at SPS. The existence of a critical density 
p c for J/ijj suppression leads to a sharp increase of t> 2 with E T , assuming that t> 2 only comes 
from J/ip dissociation in the deconfmed matter. At the RHIC energies, we found that the 
anomalous suppression already plays a role in peripheral collisions because of high density. 
In noncentral collisions there is sizable t> 2 for J/ip at pt = 3GeV. 

The experimental study of t> 2 for J/ip suppression will be complimentary to other stud- 
ies of J/ip suppression. Together with measurements, such as high pt hadrons where jet 
quenching plays an important role |24| , these studies will provide valuable information about 
the early stage of high-energy heavy-ion collisions. 
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FIG. 1. The J/ip survival probability and V2 at SPS as a function of transverse energy Er- 
nuclear absorption alone (dot-dashed line); nuclear absorption plus QGP dissociation without 
(dashed line) and with transverse expansion effects. The experimental data are from NA50 
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FIG. 2. The J ftp survival probability and v<i at RHIC as a function of number of participants 
N P . 



10 



